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Abstract 

I summarize some recent work on supersymmetric neutralinos as candidates for cold 
Dark Matter in the Universe. This includes a new scan of mSUGRA parameter space, 
with special emphasis on neutralinos annihilating predominantly through exchange of 
the light CP-even Higgs boson, and on bounds on sparticle masses. Next, prospects of 
testing models with TeV higgsino-like Dark Matter at colliders are discussed. Finally, 
I briefly comment on extensions of the mSUGRA model, and on scenarios with non- 
standard cosmology. 
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1 Introduction 



At least within the context of Einsteinian (or, indeed, Newtonian^ ) gravity, evidence for the 
existence of Dark Matter (DM) is overwhelming. This goes back to Oort and Zwicky, who 
showed in the 1930 's that there must be Dark Matter in our galaxy and in distant clusters 
of galaxies, respectively 2\. Later measurements of the total DM density, using ever larger 
systems, as well as improved determinations of the total baryon content of the Universe from 
analyses of Big Bang nucleosynthesis, implied that most of the DM must be non-baryonic. 
This was confirmed recently by analyses of data on the cosmic microwave background (CMB) 
anisotropies, in particular by the WMAP satellite, as well as other data on the large scale 
structure of the Universe, which allowed to determine many cosmological parameters with 
unprecedented precision j^J. One should be aware, however, that these determinations are 
quite indirect, and hence rely on a host of assumptions about the history of the Universe. 
In particular, one has to assume a rather simple form of the primordial power spectrum, as 
predicted by simple models of inflation. While these assumptions seem rather reasonable to 
me (see |4 j for a more contrarian view), and lead to an overall quite satisfactory (though not 
perfect) description of all existing data, it is certainly possible that these assumptions are 
not quite correct. In this talk I will therefore use a rather conservative (nominal) 3cr error 
band for the quantity of interest here, 

0.087 < n DM h 2 < 0.138. (1) 

Here Qdm is the DM mass density in units of the critical density, and h is the scaled Hubble 
constant; observations imply h 2 ~ 0.5 with ~ 10% error. The combination Qh 2 is relevant, 
since this quantity can most easily be computed theoretically. 

These calculations start from the assumption that the DM particles \ once were in full 
thermal (i.e., both kinetic and chemical) equilibrium with SM particles. This requires the 
Universe to once have been hotter than about m x /20. Since the equilibrium x number 
density is suppressed exponentially by the Boltzmann factor for T < m x , the rate T x of 
reactions that produce 1 '' or destroy DM particles would eventually have become smaller than 
the expansion rate H . At this point the DM particles are said to decouple. Under the 
assumption that the Universe evolved adiabatically after \ decoupling, one finds for today's 
DM density 

Vtv M h 2 oc . (2) 

Here (...) denotes thermal averaging, v is the relative velocity of two annihilating \ particles 
in their cms frame, and <7 e ff is the effective annihilation cross section of two \ particles into 
SM particles; in the absence of co-annihilation (see below) this is nothing but the usual 

recent attempt to construct an alternative consistent theory of gravity PP seem to me, as a particle 
physicist, far more baroque than typical particle physics models of Dark Matter, which are occasionally 
accused of being arbitrary by traditional astronomers. 

■'■Note that only SM particles with E > m x can produce DM particles; the density of such energetic SM 
particles is suppressed by the same Boltzmann factor. 
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annihilation cross section u{xx ~^ X$m), summed over all final states X$m that only contain 
SM particles. Note that the constant of proportionality in eq.(J2J) is inversely proportional 
to the Hubble parameter at decoupling. In standard cosmology, decoupling occurs when 
the Universe was radiation dominated, so that H oc g*T 2 /Mp\, where g* is the number of 
effective (relativistic) degrees of freedom and Mp\ is the Planck mass. In this case eq.(J2J) 
leads to a relic density in the range (JIJ for roughly weak-scale cross section, i.e. x should 
be a weakly interacting massive particle (WIMP). 

Under the stated reasonable assumptions - of a sufficiently hot, adiabatically evolving 
Universe containing essentially only SM particles after the DM particles decoupled - we 
have succeeded in reducing the problem of calculating the DM relic density to a problem 
of calculating cross sections. This requires choosing a particle physics model for x- Here I 
will assume that x is the lightest neutralino Xi, which is the only possible DM candidate 
in the visible sector of the Minimal Supersymmetric Standard Model 6J. From the particle 
physics point of view, this is the by far best motivated WIMP candidate. In the next Section 

1 will survey the DM-allowed regions of mSUGRA parameter space, with special focus on 
two quite different scenarios: a bino-like neutralino annihilating mostly through exchange 
of the light CP-even Higgs boson, and a heavy higgsino-like neutralino. Then I will briefly 
comment on non-minimal scenarios, before presenting my conclusions. 

2 Neutralino Dark Matter in mSUGRA 

The mSUGRA model |E] allows to describe the entire spectrum of superparticles and Higgs 
bosons by four continuous parameters and a sign, 

(m , m 1/2 , Aq, tan/?, sign(/i) ) . (3) 

Here, mo is the common (running) soft SUSY breaking scalar mass, taken at scale Mx — 
2-10 16 GeV where the SM gauge couplings meet; m 1( / 2 is the common (running) gaugino mass 
at that scale; A is the common (running) trilinear scalar soft breaking parameter, also taken 
at scale Mx', tan/? is the ratio of vacuum expectation values (vevs) of the two neutral Higgs 
fields present in the MSSM, taken at the weak scale; and \x is the supersymmetric Higgs(ino) 
mass parameter^. The parameters in (jSJ) determine the spectrum of physical masses through 
a set of renormalization group equations, as well as "threshold" loop corrections to translate 
running into pole masses. In the results presented below these are calculated using the 
Suspect program package [7j. The resulting spectrum has to pass several constraints before 
being accepted as a possible description of Nature. These include collider (chiefly LEP) limits 
on the production of superparticles and Higgs bosons, radiative b — > 37 and b — > s£ + £~ 
decays, and the DM constraint (JTJ), interpreted as described in the discussion of eq.(J2J). 
Further details can be found in [HI El - 

§The absolute value of fj, is fixed by the requirement of proper electroweak symmetry breaking. The sign 
of ji is not scale dependent. 
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2.1 Scans of parameter space 




Figure 1: The (mu2, tuq) plane in mSUGRA for tan (3 = 10 (left) and 50 (right). The values 
of the other parameters are: Aq = 0, m t = 178 GeV, and fi > 0. See the text for the meaning 
of the different shaded regions. 

Examples of scans of the mSUGRA parameter space are shown in Fig. 1. The black 
regions are excluded by theoretical constraints (in particular, by the requirement that the 
lightest superparticle, which is stable, must not be charged), as well as by the searches for 
sparticles. The pink regions are excluded by searches for neutral Higgs bosons at LEP. 
These LEP data show some (weak) evidence for the existence of an SM-like Higgs boson 
with mass near 114 GeV; the mSUGRA regions that can explain this small excess of Higgs- 
like events are shown in red. The blue region is favored by recent measurements ^U] of 
the anomalous magnetic moment of the muon, if data from hadron production at e + e~ 
colliders are used to predict the SM contribution; in this case a positive SUSY contribution 
is required at the ~ 2.5 o level. However, if instead hadronic r decay data are used to 
evaluate the SM prediction, the data are compatible with a vanishing, or even small negative, 
SUSY contribution to g^. The green region is excluded by the measured branching ratio for 
radiative b — > sj decays. As argued in [Oj, this bound is somewhat suspect; it can certainly 
be circumvented without affecting collider (or DM) phenomenology [TT]. This region is 
therefore shown behind, rather than in front of, the other regions. Finally, in the yellow 
regions the x? relic density lies in the desired range (0). The results shown in Fig. 1 broadly 
agree with those of other recent mSUGRA scans [T2"] . 

One recognizes several distinct regions with correct DM relic density. In the bulk region 
[T5) slepton masses are sufficiently small to allow a large XiXi ~ y cross section. However, 
this region is getting squeezed severely by Higgs searches at LEP. In the f co-annihilation 
region jHj the Xi~ ^i mass splitting is so small that co-annihilation processes like XiTi ~~ > T l 
have to be included in the calculation of the effective LSP annihilation cross section since 
they also reduce the total number of superparticles in the early Universe. This increases the 
effective annihilation cross section in eq.(J2J) by about one order of magnitude. These regions 
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Figure 2: The h— pole region in the (Ao,mo) plane of mSUGRA parameter space, for mt 
178 GeV, mi/2 = 140 GeV, tan/5 = 10 and \x > 0. Notation is as in Fig. 1. 



lie close to the lower black regions in Fig. I. 1 The A-pole region [13. is characterized by 
2m^o ~ rriA, so that LSP annihilation is enhanced by the near-resonant exchange of the CP- 
odd Higgs boson A in the s— channel. In mSUGRA this can only be realized at large tan/3, 
since the contribution from the b Yukawa coupling, which scales like 1/ cos/?, is essential for 
reducing through renormalization group running J7|. At large m 3> 777.1/2 one finds 
the focus point [TSj or hyperbolical branch region, where the LSP has a sizable or even 
dominant higgsino component. The location of this region depends very strongly on the 
precise value of the top mass. 

Finally, upon closer inspection you will see small DM-allowed regions at tuq ~ 2 to 3 
TeV and close to its lower bound (which essentially comes from chargino searches at 
LEP). In these regions of parameter space the Xi annihilation cross section is enhanced by 
the near-resonant exchange of the light CP-even Higgs boson h. This region has therefore 
been called the h—pole region in It featured prominently in pre-LEP2 mSUGRA scans 
(see e.g. [T3J 1201), but then was squeezed almost out of existence by Higgs and chargino 
searches. It has recently been resurrected [21] by improved calculations of m^, as well as by 
the somewhat increased central value of the top mass (which, however, shows a very recent 
tendency of falling again, if preliminary data j22] are to be trusted). For m t = 178 GeV, the 
h— pole region requires 

mi/2 e [130 GeV, 145 GeV] (h - pole region). (4) 
conceptually similar stop co-annihilation region ^B] is difficult to realize in mSUGRA. 
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However, as shown in Fig. 2, it covers quite a substantial area in the (Aq, m ) plane for 
m 1/2 = 140 GeV. 

2.2 Mass bounds 

Bounds on physical (pole) masses might be a more meaningful way to show the possibilities of 
mSUGRA than the ubiquitous plots of allowed regions in the space of basic input parameters 
(JHJ). For one thing, these plots always fix the values of other free parameters. Usually Aq = 
is chosen, but nonvanishing values of Aq can have quite significant impact on allowed spectra 
[HI 123 • Similarly, the current error on the top quark mass is still significant. Varying m t 
in its 90% c.l. interval (according to [24. ) moves around the lower bounds on m or 
from Higgs searches by hundreds of GeV; the shift of the focus point / hyperbolical branch 
region is even bigger. Finally, the parameter tan (3 is still almost completely unconstrained, 
apart from the lower bound tan/3 > 1.5 that follows from Higgs searches in mSUGRA (for 
stop masses below 2 TeV). As well known ^7j, the sparticle spectrum for models with large 
tan j3 looks quite different from that for low or moderate values of tan j3, with quite drastic 
impact on collider physics |2H] . 

One obtains the least biased view of the allowed ranges of masses by simply scanning 
over the entire parameter set that is consistent with a given set of constraints." Not all 
constraints should be treated on an equal footing. Lower bounds on masses (or cross sections 
or branching ratios) from accelerator-based experiments are most robust, since both beam 
and detector are well controlled by the experimenter (I hope). Bounds on masses and 
tree-level cross sections usually also don't have many theoretical ambiguities. Bounds on 
processes that only occur at quantum level, like b — > decays or the anomalous magnetic 
moment of the muon, are already somewhat less robust, since here in general several diagrams 
contribute, which opens the possibility of cancellations. More generally, one can often evade 
those bounds by relatively minor modifications of the model. In case of — 2, there's 
an additional ambiguity due to the ~ 2a discrepancy between SM 'predictions' based on 
different data sets, as mentioned in the discussion of Fig. 1. 

Including the DM constraint (JIJ is fine in principle, since the existence of DM really has 
been established beyond reasonable doubt in the framework of mSUGRA or similar models.** 
However, translating this into a bound on model parameters, or on physical masses, is far 
from trivial. The usual procedure [T7| |2UJ |H1 C2J is to take eq.(J2J) at face value. However, as 
discussed in the Introduction, this means that one accepts a set of assumptions about the 
early Universe that are presently untested, and indeed are difficult to test except through 
Dark Matter physics. Nevertheless, since these assumptions are indeed reasonable, it makes 
sense to try and see where they lead us. 

"i do not see much value in attempts to assign probabilities to certain regions of parameter space. 
These parameters cannot be interpreted as randomly distributed quantities, unless one subscribes to "land- 
scape" and "multiverse" ideas, which to my mind come dangerously close to pseudoscience. If these param- 
eters are indeed fixed, it makes little sense to assign probability distributions to them. 

**If one entertains the idea of modifying gravity 1 rather than introducing Dark Matter, supergravity 
would certainly also have to be modified drastically. 
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It was hoped originally that the upper bound on the DM relic density (the so-called 
'overclosure' constraint) would allow to establish reliable, useful upper bounds on sparticle 
masses. Eq. (J2J) shows that the relic density goes like the inverse of the (effective) Xi annihi- 
lation cross section, which in turn (through dimensional arguments, or by unitariy goes 
like the inverse square of the relevant mass scale. Indeed, unitarity does allow to establish 
an upper bound on the mass of any WIMP; however, this bound exceeds 100 TeV [27], and 
is thus not particularly useful, since we lack the means to build colliders that could cover 
this kind of mass range. In the context of mSUGRA, it became clear quite early on [13] that 
very, even "unnaturally", large masses can be compatible with the DM constraint (0) even 
in standard cosmology. 
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Table 1: Sparticle mass bounds in mSUGRA. These have been obtained by scanning over the 
entire allowed parameter space, defined by m t G [171 GeV, 185 GeV], (m^+mj )/2 < 2 TeV, 
the lower bounds on sparticle and Higgs masses from collider experiments, the constraint on 
b — *■ sj implemented as in jH], simple 'CCB' constraints j2Hj, and a conservative interpretation 
of the constraint from — 2 (essentially the overlap of the 2a regions using r decay and e + e~ 
collider data [9J. Set II adds the DM constraint (JTJ), implemented using (J2J), to the above 
set of constraints. Set III is like Set II, except that the scanned region as been artifically 
limited to the h— pole region, where m^o < Only lower bounds are listed for Sets I 

and II, while for Set III the allowed range is given; a dash (-) means that the upper bound 
is directly set by the arbitrary upper bound on the average stop mass. 

Nevertheless it is clear from Figs. 1 and 2 that this constraint does exclude very large 
chunks of otherwise allowed parameter space. One might therefore think that it would 
at least affect the lower bounds on sparticle masses significantly. The Table below shows 
that this is not really the case. This table lists lower bounds on the masses of some new 
(s)particles in mSUGRA, first without (set I) and then with (set II) including the DM 
constraint. In particular, the lower bounds on many new (s)particles simply coincide with 
the bounds established by collider experiments. This is true for the lighter chargino, stau and 
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scalar Higgs states, and essentially also holds true for the lighter stop. The bounds on the 
masses of the gluino and third neutralino are essentially the same as that in a more general 
MSSM, as long as gaugino mass unification is maintained. Clearly the DM constraint still 
allows some new (s)particles to be quite light! I should emphasize, however, that usually the 
lower bounds in the Table cannot be saturated simultaneously. E.g. the lower bounds on 
electroweak gaugino masses are saturated at > 1 TeV, where first generations squarks 
and sleptons are quite heavy. Nevertheless, the possibility of light sparticles even in this 
simplest of all potentially realistic SUSY models that allow WIMP Dark Matter should be 
quite encouraging to experimenters! 

Set III shows these bounds (including the DM constraint) when one confines oneself 
to the h— pole region discussed at the end of the previous subsection. In this case there 
are significant upper bounds on the masses of all gauginos. The reason is that one needs 
2m^o ~ m h < 120 GeV here. Since m^o ~ m~± ~ M 2 — 2M 1 and m~ g ~ M 3 ~ 6M X in 
mSUGRA (due to the assumed universality of gaugino masses at the GUT scale), one finds 
quite stringent upper bounds on the masses of the lighter chargino and gluino. These bounds 
can be tested at the Tevatron, if it manages to collect several fb _1 of data; testing the gluino 
mass bound at the LHC is essentially trivial. On the negative side, we see that most scalars 
have to be quite heavy in this region. This is required to satisfy the h search limit from LEP 
with small m^. 

2.3 Heavy higgsino Dark Matter 

In the last subsection I imposed the naturalness constraint that the average stop mass should 
not exceed 2 TeV. If this constraint is relaxed, the deep hyperbolical branch region becomes 
allowed. For m ~ 2m 1 / 2 > 6 TeV (for m t = 178 GeV) one finds [2H] DM-allowed points 
where the LSP is an almost pure higgsino with mass near 1 TeV [3U]. This region is impossible 
to probe even at the LHC, since squark and gluino masses lie near 10 TeV or even higher; 
in fact, mo and could even be exponentially larger than the weak scale! A decisive 
test will therefore have to make do with the three relatively light sparticles, i.e., the three 
higgsino-like states Xii Xi an d X.2- Currently the best bet for producing these sparticles in 
detectable quantities seems to be the futuristic e + e~ collider CLIC jSI], which may operate 
at center-of-mass energy y/s = 3 TeV. However, since the mass splittings between these 
higgsino states amount to at most a few GeV, detecting them will not be straightforward. 
The reason is that one expects large fluxes of beamstrahlung photons at future e + e~ colliders, 
which tend to get worse with increasing beam energy and increasing luminosity. (Note that 
the luminosity should exceed that of LEP by a factor of order sclic / s lep ~ 200 to produce a 
useful number of e + e~ annihilation events.) This flux might be so large that each annihilation 
event is accompanied by an "underlying event" consisting of one or more soft 77 — > hadrons 
reactions [3*2*] . 

This can be overcome by demanding that the higgsinos are produced together with a 
hard photon, that is emitted at large angle to the beam. A background event would then 
have to contain a hard e + e~ — > uPj annihilation plus an underlying event that happens to 
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look like the soft x° or xt decay products. A recent calculation [22 EH] indicates that this 
background (and similar ones occurring at higher orders in electroweak couplings) can be 
overcome, especially if longitudinally polarized beams are available. 

3 Nonminimal scenarios 

The small size of the DM-allowed parameter space in Fig. 1 is occasionally taken as moti- 
vation to consider extensions of the mSUGRA model. Actually, a small allowed parameter 
space per so should not be considered to be problematic. New measurements, like the con- 
straint are bound to constrain the parameter space of any theory. In fact, ideally one 
wants to (over-) constrain the parameter space, so that it collapses to (almost) a point, as in 
QED (considered here as interaction of electrons with photons, i.e. a 2 parameter theory). 

Slightly more worrisome is perhaps the observation that almost everywhere in the DM- 
allowed mSUGRA parameter space some particular relation between parameters is required 
to hold (in the pole or co-annihilation regions), and/or one is close to the edge of the region 
excluded by consistency considerations: in the f co-annihilation region one is close to the 
region where the T\ would be LSP, leading to a forbidden charged stable particle, whereas 
the focus point / hyperbolical branch region is close to the edge of parameter space where 
consistent electroweak symmetry breaking no longer works. One might therefore ask whether 
more 'generic' DM-allowed regions of parameter space exist in (slight) generalizations of 
mSUGRA. 

The answer, I believe, is No. This is illustrated by Fig. 3, which shows the scaled DM 
relic density as function of one parameter, with all other parameters held fixed. Starting 
point is an mSUGRA scenario with m t = 178 GeV, m = 800 GeV, m 1 / 2 = 500 GeV, A Q = 0, 
tan/3 = 10 and /i > 0. On the red curve, which is actually still an mSUGRA curve, A has 
been increased. This reduces the mass of the lighter stop eigenstate ti, eventually ti — Xi 
co-annihilation [16. sets in, reducing the relic density. Alas, f^o/i 2 quickly undershoots the 
desired range, i.e. one gets the right relic density only for a very narrow range of values of 

The situation is similar on the blue line, where I have increased both Higgs soft breaking 
masses at the GUT scale relative to the sfermion masses [33] . This allows to reduce the value 
of at the weak scale, making the LSP more higgsino-like. However, again the relic density 
quickly falls below the desired range; as we saw in the last Subsection of the previous Section, 
a higgsino-like LSP would need a mass near 1 TeV to make a good thermal CDM candidate; 
this is well above the LSP masses considered in Fig. 3. The DM constraint therefore again 
singles out a very narrow strip of the (extended) parameter space. Note also that in this, 
as well as the previous, case the DM-allowed region is again quite close to the upper bound 
on the parameter that is being varied (due to the experimental bound on and on m ~± , 
respectively) . 

The dark green curve is for a model where gaugino masses are not universal at the GUT 
scale [35]. Specifically, I have kept Mi = M 3 , but reduced M 2 , making the LSP more 
wino-like. Unfortunately the curve yet again falls quickly through the desired range of relic 
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Figure 3: The scaled relic density in slight generalizations of mSUGRA. The desired range 
is indicated by the light green band. Demanding that the relic density falls in this band 
singles out a narrow sliver of the otherwise allowed range of the parameter P that is being 
varied. See the text for further details. 



density, i.e. the DM constraint is again satisfied only in a narrow sliver of parameter space; 
in this case the allowed region is not close to the edge, though. 

On can also consider generalizing the cosmological framework. All curves in Fig. 3 still 
use the standard assumptions (sufficiently high temperature, only SM particles contribute 
to the Hubble expansion, no entropy production after decoupling of the LSP). For example, 
a quintessence field could significantly alter the Hubble parameter during decoupling. As 
far as I know, only scenarios with increased H have been analyzed so far [HE]. In this case 
one can increase the relic density by a large factor, e.g. allowing higgsino Dark Matter with 
"natural" mass of a few hundred GeV. However, one needs to tune the model to have a large 
change of H during LSP decoupling and in the present period (to account for the observed 
accelerated expansion of the Universe), while keeping H very close to its standard value 
during Big Bang nucleosynthesis. 

Another way to increase the LSP density is to assume that there were non-thermal 
production mechanisms [37]. Finally, if the lightest neutralino is long-lived but eventually 
decays into a very weakly coupled superparticle, e.g. an axino j3S] or gravitino jHU], one 
automatically becomes more sensitive to details of early universe cosmology. The simplest 
assumption is that the decay of the neutralinos (or whatever else is the lightest sparticle 
in the visible sector) is the by far dominant contribution to the relic density of this super- 
weakly interacting sparticle [HI]- However, in general there can be significant production 



10 



of gravitinos or axinos also from the thermal bath [II]. Finally, there might again be ad- 
ditional non-thermal production mechanisms [12] • Scenarios with very weakly interacting 
dark matter particles therefore allow to make any point in the parameter space of mSUGRA 
(or similar models) DM-compatible. If the lightest neutralino remains the lightest visible 
sparticle, these models often have exactly the same collider phenomenology as if x? indeed 
was the LSP. 

4 Summary and conclusions 

The constraint (0) greatly constrains the parameter space of any particle physics model 
of Dark Matter. The mSUGRA model still works fine, although some relation between 
parameters seems to be required. One (resurrected) possibility is 2m^o ~ rrih (the h— pole 
region). More generally, in mSUGRA the DM constraint does not increase lower bounds on 
sparticle and Higgs masses significantly. Very (almost arbitrarily) large gaugino and sfermion 
masses are allowed if the LSP is higgsino-like; it might nevertheless be possible to test this 
region of parameter space at CLIC. Generalizations of mSUGRA allow to reduce the DM 
relic density, but generally do not increase the fraction of DM-allowed parameter space. 
Finally, if the LSP resides in the hidden sector, pretty much any visible sector spectrum 
can be made DM-safe, at the cost of predictivity and testability. I conclude that the Dark 
Matter constraint is a great tool for model builders, but should probably not be taken too 
seriously by collider physicists. 
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